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Operator fractional Brownian motions (OFBMs) are (i) Gaussian, (ii) operator self-similar and 
(iii) stationary increment processes. They are the natural multivariate generalizations of the 
well-studied fractional Brownian motions. Because of the possible lack of time-reversibility, the 
defining properties (i)-(iii) do not, in general, characterize the covariance structure of OFBMs. 
To circumvent this problem, the class of OFBMs is characterized here by means of their inte- 
gral representations in the spectral and time domains. For the spectral domain representations, 
this involves showing how the operator self-similarity shapes the spectral density in the general 
representation of stationary increment processes. The time domain representations are derived 
by using primary matrix functions and taking the Fourier transforms of the deterministic spec- 
tral domain kernels. Necessary and sufficient conditions for OFBMs to be time-reversible are 
established in terms of their spectral and time domain representations. It is also shown that the 
spectral density of the stationary increments of an OFBM has a rigid structure, here called the 
dichotomy principle. The notion of operator Brownian motions is also explored. 

Keywords: dichotomy principle; integral representations; long-range dependence; multivariate 
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1. Introduction 

Fractional Brownian motion (FBM), denoted Bh ~ {BH{t)}teM with H e (0,1), is a 
stochastic process characterized by the following three properties: 

(i) Gaussianity; 

(ii) self-similarity with parameter H; 

(iii) stationarity of the increments. 

By self-similarity, it is meant that the law of Bh scales as 

{BHict)}tm^{c"BH{t)}tm, c>0, (1.1) 
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where = denotes equality of finite-dimensional distributions. By stationarity of the in- 
crements, it is meant that the process 



has the same distribution for any time-shift /i G R. It may be shown that these three 
properties actually characterize FBM in the sense that it is the unique (up to a constant) 
such process for a given H £ (0, 1). FBM plays an important role in both theory and 
applications, especially in connection with long-range dependence [11, 12]. 

We arc interested here in the multivariate counterparts of FBM, called operator 
fractional Brownian motions (OFBMs). In the multivariate context, an OFBM Bh = 
{Bi,H, ■ ■ ■,Bn,H)* = {{Bi.H{t),- ■ .,Bn,H{t))* G R'N^ G M} is a collection of random vec- 
tors, where the symbol * denotes transposition. It is also Gaussian and has stationary 
increments. Moreover, as is standard for the multivariate context, in this paper, we as- 
sume that OFBMs are proper, that is, for each t, the distribution of Bnit) is not contained 
in a proper subspacc of R". However, self-similarity is now replaced by 

(ii') operator self-similarity. 

A proper multivariate process Bh is called (strictly) operator self-similar (o.s.s.) if it 
is continuous in law for all t and the expression (1.1) holds for some matrix H . Here, the 
expression is defined by means of the convergent series 



Operator self-similarity extends the usual notion of self-similarity and was first studied 
thoroughly in [19, 22]; see also Section 11 in [28] and Chapter 9 in [12]. The theory of 
operator self-similarity runs somewhat parallel to that of operator stable measures (see 
[20, 28]) and is also related to that of operator scaling random fields (see, e.g., [3]). 

OFBMs arc of interest in several areas and for reasons similar to those in the univari- 
ate case. For example, OFBMs arise and are used in the context of multivariate time 
scries and long-range dependence (sec, e.g., [4-6, 9, 26, 30]). Another context is that 
of queueing systems, where reflected OFBMs model the size of multiple queues in par- 
ticular classes of queueing models and are studied in problems related to, for example, 
large deviations (see [7, 21, 24, 25]). Partly motivated by this interest in OFBMs, several 
authors consider constructions and properties of OFBMs. Maejima and Mason [23], in 
particular, construct examples of OFBMs through time domain integral representations. 
Mason and Xiao [27] study sample path properties of OFBMs. Bahadoran, Benassi and 
Dgbicki [1] provide wavelet decompositions of OFBMs, study their sample path proper- 
ties and consider questions of identification. Becker-Kern and Pap [2] consider estimation 
of the real spectrum of the self-similarity exponent. A number of other works on opera- 
tor self-similarity are naturally related to OFBMs; see, for example, [28], Section 11, and 
references therein. 

To the reader less familiar with OFBMs, we should note that the multivariate case 
is quite different from the univariate case. For example, consider an OFBM Bh whose 



{BH{t + h)-BH{h)}tm 
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exponent H has characteristic roots hk with positive real parts. By using operator self- 
similarity and stationarity of increments, one can argue, as in the univariate case, that 

EBH{t)BH{sy + EBH{s)BH{t)* 

= EBH{t)BH{ty + EBh{s)Bh{s)* - E{BH(t) - BH{s))(BH{t) - Bh{s))* (1.2) 

= \t\"T{i,i)\t\"' + \s\"r{i,i)\s\"' -\t-s\"rii,i)\t-s\"', 

where r(l, 1) = EB h{1)B h{1)* and the symbol * denotes the adjoint operator. However, 
in contrast with the univariate case, it is not generally true that 

EBH{t)BH{sy = EBH{s)BH{t)* (1.3) 

and hence the OFBM is not characterized by H and a matrix r(l,l). Another impor- 
tant difference is that the self-similarity exponent of an operator self-similar process is 
generally not unique. The latter fact has been well known since the fundamental work 
of Hudson and Mason [19]. We briefly recall it, together with some related results, in 
Section 2.2. 

In this work, we address several new and, in our view, important questions about 
OFBMs. In view of (1.2) and (1.3), since the covariance structure of an OFBM cannot 
be determined in general, we pursue the characterization of OFBMs in terms of their 
integral representations (Section 3). In the spectral domain, under the mild and natural 
assumption that the characteristic roots of H satisfy 

0<Re(/ifc)<l, fc = l,...,n, (1.4) 

we show that an OFBM admits the integral representation 

{x-^"-™A + xZ}"~^"^^'^A)B{Ax), (1.5) 

where A is a matrix with complex-valued entries, A denotes its complex conjugate, 
a;+ = max{a;,0}, x_ = max{— a;,0} and B{Ax) is a suitable multivariate complex-valued 
Gaussian measure. In the time domain and when, in addition to (1.4), we have 

Re(/ife)^i, fc = l,...,n, (1.6) 
the OFBM admits the integral representation 

[ ((1+ \H-{l/2)I , ^^r-(l/2)7^ , ,r , ■,H-[l/2)I , ^^r-(l/2)7^ , r X 

/ {{[t~u)_^ - (-w)^ )M+ + ((t-ii)_ ~[~u)_ )M^)B(du), 

Jr 

(1-7) 

where Af-j., A/_ are matrices with real- valued entries and B{du) is a suitable multivariate 
real- valued Gaussian measure. The representation (1.7) is obtained from (1.5) by taking 
the Fourier transform of the deterministic kernel in (1.5). We shall provide rigorous 
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arguments for this step by using primary matrix functions. (Even in the univariate case, 
very often this step is unjustifiably taken as more or less evident.) On a related note, 
but from a different angle, the representations (1.5) and (1.7) always define Gaussian 
processes with stationary increments that satisfy (1.1) for a matrix H. We shall provide 
sufficient condition for these processes to be proper (see Section 4) and, hence, to be 
OFBMs. 

Subclasses of the representations (1.5) and (1.7) were considered in the works refer- 
enced above. Maejima and Mason [23] consider OFBMs given by the representation (1.7) 
with = A/_ = /. Mason and Xiao [27] take (1.5) with I. Bahadoran et al. [1] con- 
sider (1.5) with A having full rank and real-valued entries. (Such OFBMs, for example, 
are necessarily time-reversible; see Theorem 5.1 and also Remark 3.1.) We would again 
like to emphasize that, in contrast with these works, the representations (1.5) and (1.7) 
characterize all OFBMs (under the mild and natural conditions (1.4) and (1.6)). 

In particular, the representations (1.5) and (1.7) provide a natural framework for the 
study of many properties of OFBMs. In this paper, we provide conditions in terms of A 
in (1.5) (or M+,M_ in (1.7)) for OFBMs to be time-reversible (see Section 5). Time- 
reversibility is shown to be equivalent to the condition (1.3) and hence, in view of (1.2), 
corresponds to the situation where the covariance structure of the OFBM is given by 

EBH{t)BH{sy = i(irr(i, \)\tr + 1)1^1^' - \t - s|^r(i, i)|^ - sf). (i.s) 

Another interesting and little-explored direction of study of OFBMs is their uniqueness 
(identification). This encompasses the characterization of the different parameterizations 
for any given OFBM and, in particular, of the aforementioned non-uniqueness of the 
self-similarity exponents. Uniqueness questions in the context of OFBMs are the focus 
of Didier and Pipiras [8], where they are explored starting with the representation (1.5), 
and will be largely absent from this paper. 

Furthermore, in this paper, we also discuss some additional properties of OFBMs 
which are of independent interest. First, we prove that OFBMs have a rigid dependence 
structure among components, which we call the dichotomy principle (Section 6). More 
precisely, under long-range dependence (in the sense considered in Section 6), we show 
that the components of the increments of an OFBM are either independent or long- 
range dependent, that is, they cannot be short-range dependent in a non-trivial way. 
Since, in the univariate case, the increments of FBM are often considered representative 
of all long-range dependent series, this result raises the question of whether OFBMs 
are flexible enough to capture multivariate long-range dependence structures. Second, 
we also discuss the notion of operator Brownian motions (OBMs) and related questions 
(Section 7). OBMs are defined as having independent increments and are known to admit 
H = (1/2)/ as an exponent. We also show, in particular, that an OFBM with H = (1/2)/ 
does not necessarily have independent increments and hence is not necessarily an OBM. 
(In contrast, in the univariate case, H = 1/2 necessarily implies Brownian motion.) 

In summary, the structure of the paper is as follows. In Section 2, wc provide the 
necessary background for the paper and some definitions. In Section 3, wc construct in- 
tegral representations for OFBMs in the spectral and time domains. Section 4 furnishes 
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sufficient conditions for properness. Section 5 is dedicated to time-reversibility. The di- 
chotomy principle is established in Section 6 and the properties of OBMs are studied in 
Section 7. Appendices A-D contain several important technical results used throughout 
the paper, as well as some proofs. 

2. Preliminaries 

We begin by introducing some notation and considering some preliminaries on the expo- 
nential map and operator self-similarity that are used throughout the paper. 

2.1. Some notation 

In this paper, the notation and terminology for finite-dimensional operator theory will 
be preferred over their matrix analogs. However, whenever convenient, the latter will be 
used. 

All with respect to the field R, M{n) or M{n,R) is the vector space of all n x n 
operators (endomorphisms) , GL(n) or GL(n, M) is the general linear group (invertible 
operators, or automorphisms), 0(n) is the orthogonal group of operators O such that 
OO* = / = 0*0 (i.e., the adjoint operator is the inverse), SO{n) C 0{n) is the special 
orthogonal group of operators (rotations) with determinant equal to 1 and so{n) is the 
vector space of skew-symmetric operators (i.e., A* = — A). 

The notation will indicate a change to the field C. For instance, M(ri,C) is the vector 
space of complex endomorphisms. Whenever it is said that A G M{n) has a complex 
eigenvalue or eigenspace, one is considering the operator embedding M(n) ^ M(n,C). 
The notation A indicates the operator whose matrix representation is entrywise equal to 
the complex conjugates of those of A. We will say that two endomorphisms A,B£ M{n) 
are conjugate (or similar) when there exists P G GL(n) such that A — PBP~^. In this 
case, P is called a conjugacy. The expression diag(Ai, . . . , A„) denotes the operator whose 
matrix expression has the values Ai, . . . , A„ on the diagonal and zeros elsewhere. The ex- 
pression tr(A) denotes the trace of an operator A £ M{n, C). We write / € L^(]R, M(n, C)) 
for a matrix- valued function / when tr{Jg /(m)*/(m) du} < +oo. 

Throughout the paper, we set 

D = H- {1/2)1 (2.1) 
for an operator exponent H . The characteristic roots of H and D arc denoted 

hk,dk, (2.2) 

respectively. Here, 

k^l,...,N or n, (2.3) 

where A'^ < n is the number of different characteristic roots of H . 

For notational simplicity when constructing the spectral and time domain filters, we 
will adopt the convention that = £ M{n, M) when z = 0. 
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2.2. Operator self-similar processes 

Operator self-similar (o.s.s.) processes were defined in Section 1. Any matrix H for which 
(1.1) holds is called an exponent of the o.s.s. process X. The set of all such H for X 
is denoted by £{X), which, in general, contains more than one exponent. The non- 
uniqueness of the exponent H depends on the symmetry group Gi of X, which is defined 
as follows. 

Definition 2.1. The symmetry group of an o.s.s. process X is the set Gi of matrices 
A G GL{n) such that 

{X{t)}te^^{AX{t)}tm- (2.4) 

It turns out that the symmetry group Gi is always compact, which implies that there 
exists a closed subgroup Oq of 0{n) such that Gi = WOqW^^, where is a positive 
definite matrix (see, e.g., [20], Corollary 2.4.2, page 61). A process X that has maximal 
symmetry, that is, such that Gi = WO{n)W^^ , is called elliptically symmetric. 

Let G be a closed (sub)group of operators. The tangent space T{G) of G is the set of 
A e M{n) such that 

G — / 

A = lim — — for some {G„} C G and some < (i„ — t' 0. 

n— >oo dn 

In this sense, T{G) is, in fact, a linearization of G in a neighborhood of /. Hudson and 
Mason [19], Theorem 2, shows that for any given o.s.s. process X with exponent H , the 
set of exponents £{X) has the form £{X) H + T{Gi), where T{Gi) = WCqW'^ for 
the positive definite conjugacy matrix W associated with Gi and some subspace Cq of 
so{n). Consequently, X has a unique exponent if and only if Gi is finite. 

3. Integral representations of OFBMs 

Representations of OFBMs in the spectral domain are derived in Section 3.1. The corre- 
sponding representations in the time domain are given in Section 3.2. The derivation of 
these representations is quite different from that in the univariate case. In the latter case, 
it is enough to "guess" the form of the spectral representation and to verify that it gives 
self-similarity and stationarity of the increments (and hence, immediately, FBM). In the 
multivariate case, these representations actually have to be derived from the properties 
of OFBMs, without any guessing involved. 

3.1. Spectral domain representations 

In Theorem 3.1, we establish integral representations of OFBMs in the spectral domain. 
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Theorem 3.1. Let H Cz M{n,]&.) with characteristic roots hk satisfying 

0<Re(/ifc)<l, fc = l,...,n. (3.1) 

Let {i?^f(i)}tgR he an OFBM with exponent H. Then {-Bh (i)}teR admits the integral 
representation 

{BHmtm - ( / ^——^{x^^'A + xZ''A)B{dx)] (3.2) 

for some ACz M{n,C). Here, D is as in (2.1), 

B{x):^Bi{x)+iB2{x) (3.3) 

denotes a complex-valued multivariate Brownian motion such that Bi{~x) = Bi{x) and 
B2{—x) = —B2{x), Bi and B2 are independent and the induced random measure B(dx) 
satisfies EB{dx)B{dx)* ^dx. 

Proof. For notational simplicity, set X = Bh- Since X has stationary increments, we 
have 



X{t)-X{s)^ / : Y{dx), (3.4) 

where Y{dx) is an orthogonal- increment random measure in C". The relation (3.4) can 
be proven following the approach for the univariate case found in [10], page 550, under 
the assumption that E\X{t + h) — X{t)\'^ — > as /i — > 0, that is, X is L^-continuous at 
every t (see also [35], page 409, and [34], Theorem 7). The latter assumption is satisfied 
in our context because of the following. Property 2.1 in [23] states that, for an o.s.s. 
process Z with exponent H , if inf{Re(/ifc); fc = 1, . . . , n} > 0, then Z{Q) = a.s. Thus, in 
view of (3.1), X[Q) = a.s. So, by stationarity of the increments and continuity in law, 

X{t + h)- X{t)^X{h)^ X{Q)=0, h^O. (3.5) 
Therefore, by relation (3.4) and again by Property 2.1 in [23], 



X{t)^ / — : Y(dx). (3.6) 

Js. 

Let Fx{dx) — EY {dx)Y (dx)* be the multivariate spectral distribution of Y{dx). The 
remainder of the proof involves three steps: 

(i) showing the existence of a spectral density function fx{x) = Fx{dx) / dx; 

(ii) decorrelating the measure Y{dx) componentwise by finding a filter based upon 
the spectral density function; 

(iii) developing the form of the filter. 
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Step (i). Since X is o.s.s. with exponent H, 

X{ct)^c" [ ^—T^Yidx) (3.7) 

for c > 0. On the other hand, through a change of variables x = c~^v, 

X{ct)= [ ~^ cY{c~^dv). (3.8) 

The relations (3.7) and (3.8) provide two spectral representations for the process 
{X(rf)}tgR. As a consequence of the uniqueness of the spectral distribution func- 
tion of the stationary process {X{t) — X{t — l)}tgR and of the fact that | ° . ~^ p > 0, 
X e K\{27tfc, e Z}, we obtain that 

c^Fx{c-^dx)=c"Fx{dx)c"\ c> 0. 

Equivalently, by a simple change of variables, Fx{cdx) = c^^^ Fx{dx)c'^^~^^ . Thus, for 
OO, 

/ Fx(cdx)=Fx(0,c]=c^-^Fx(0,l]c(^-^)*, (3.9) 
J(04] 

^x(cdx) =Fx(-c,0] =c^"^Fx(-l,0]c(-^--^'*. (3.10) 

(-1.0] 

By the explicit formula for c^~^ in Appendix D. each individual entry _F\'(0,c]ij, i,j = 
1, . . . ,ri,, in the expression on the right-hand side of (3.9) is either a linear combination 
(with complex weights) of terms of the form 

(log(c))' i_,, Jlog(c))" 
/! to! 

or is identically zero for c > 0. Thus, Fx{c) is differentiable in c over (0,00) since 
Fx(0,c]y = Fx{c)ij-Fx{0)ij. The differentiability of Fx on (-00, 0) follows from (3.10) 
and an analogous argument. 

To finish the proof of the absolute continuity of Fx, it suffices to show that Fx is 
continuous at zero. Note that 

Fx(-c,c] =c'^"-^Fx(-l,l]c('^"-^'' ^0 

as c — 0+. The limit holds because ||c^^^|| — ^ as c — 0"*", where || • || is the matrix 
norm, which, in turn, follows from [23], Proposition 2.1(ii), under the assumption that 
Re{hk) < 1, fc= l,...,n. 

Step (ii). Denote the spectral density of X hy fx - Since I''— ^ — \^fx{x) is the spectral 
density of the stationary process {X{t) — X{t — l)}t£R, fx{x) is a positive semidefinite 
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Hermitian symmetric matrix dx-a.e. ([15], Theorem 1, page 34). The spectral theorem 
yields a (unique) positive semidefinite square root a{x) of fx{x). Let B{x) be a complex- 
valued multivariate Brownian motion, as in the statement of the theorem. X can then 
also be represented as 

X{t)^ I '—T^Z{x)B{dx) (3.12) 



because 

E{Z{x)B{Ax)B{Ax)*a{x)*)^a{xf Ax = fx{x) dx = Fx (da;) 

and the processes on both sides of (3.12) are Gaussian and real-valued. 

Step (iii). By using operator self-similarity and arguing as in step (i), the relation 
(3.12) implies that, for every c > 0, 

-D^( ^ \^( X \ _£3* 



a{x)a{x)* ~ c ay— ja\^— j c dx-a.e. (3.13) 

By Fubini's theorem, the relation (3.13) also holds dxdc-a.e. 
Consider x > 0. A change of variables leads to 

a{x)a{x)* ~ x~^v^a{v)a{v)*v^ x~^ dxdv-a.e. 
Thus, one can choose u+ > such that 

a{x)a{xy =x-°v^a{v+)a{v+yvfx-°' dx-a.e. (3.14) 

This means, in particular, that if we set 

S+(x) = x~^v^a{vj^-) 

for dx-a.e. a; > 0, then a^{x)a+{x)* = fx{x) on the same domain. 

By again considering the stationary process {X{t) — X{t— l)}teR and applying [15], 
Theorem 3, page 41, one can show that fx is a Hermitian function. Thus, 



a{—x)a{~x)*=fx{—x) = fx{x)=x v^a{v+)a{v+)*v^ x dx-a.e. 
Hence, for x < 0, we can set 



a_(x) = {—x)^^ v^a{v+) 

and, for x G M, we have 



5;(x) = x_|_ v^a{vj^) + x_ v^a{vj^) dx-a.c., 

where a(x)a(x)* = fx{x) dx-a.e. Therefore, we can use a in place of a in the spectral 
representation of X, which establishes relation (3.2). □ 
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Remark 3.1. The invertibility of A in relation (3.2) is not a requirement for the process 
to be proper (compare with [1], page 9). In the Gaussian case, properness is equivalent 
to EX{t)X{t)* being a full rank matrix for all t^Q. 

A simple example would be that of a bivariate OFBM whose spectral representation 
has matrix parameters D = dl, < d < 1/2, and A set to the (unique) non-negative 
square root of 

which is rank-deficient. Let 



/>oo 


cit^ _ 1 


Jo 


ix 



9{t) 

which is strictly positive for all t^Q. In this case, 

2 



-2d 



da;, 



Ex{t)x{ty 



1 



9{t) 



IX 

2 _ 

i + i 



|xr^''(A^l{,>o}+A2l{,<o})d.T 
* + 2g{t)I. 



Theorem 3.1 shows that an OFBM is characterized by a (potentially non- unique) 
o.s.s. exponent H and a matrix A. For the sake of simplicity, we will continue to use the 
notation Bh instead of the (more correct) notation Bh,a- 

Remark 3.2. As a consequence of Macjima and Mason [23], Corollary 2.1, the charac- 
teristic roots hk of the exponent H of an OFBM must satisfy Rc(/ifc) <1, A: = l,...,n. 
However, the extension of the definition of OFBMs to the case of H with at least one 
characteristic root hk satisfying Re(/ifc) ~ 1 can be subtle. In Proposition C.l, it is shown 
that there does not exist an OFBM with exponent 




whose characteristic roots are = /12 = 1. (More precisely, it is shown that a Gaus- 
sian, H-o.s.s. process X = {Xi,X2)* with stationary increments is necessarily such that 
Xi{t) = and X2{t) = tY a.s. for a Gaussian variable Y and hence that it cannot be 
proper.) 



3.2. Time domain representations 

Our next goal is to provide integral representations of OFBMs in the time domain, 
which is done in Theorem 3.2. The key technical step in the proof is the calculation of 
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the (entrywise) Fourier transform of the kernels 

{t - u)^ - i-u)^ = cxp(log(t - u)±D) - cxpi\og{-u)±D), (3.15) 

which are the muhivariate analogs of the corresponding univariate FBM time domain 
kernels. It is natural and convenient to carry out this step in the framework of the so- 
called primary matrix functions. The latter allows one to naturally define matrix analogs 
f{D), D £ M (n, R), of univariate functions ,f{d), d G M, and to say when two such matrix- 
valued functions are equal based on their univariate counterparts. 

For the reader's convenience, we recall the definition of primary matrix functions (more 
details and properties can be found in [17], Sections 6.1 and 6.2). Let A e Af(n, C) with 
minimal polynomial 

qA{z)^{z~XiY'---{z~XNr-, (3.16) 

where Xi,. . . ,Xn are pairwisc distinct and > 1 for k ^ 1, . . . , N , N <n. We denote by 
A = PJP~^ the Jordan decomposition of A, where J is in Jordan canonical form with 
the Jordan blocks Jx^ , . . . , on the diagonal. 

Let J7 C C be an open set. Given a function h:U ^ C and some A G M{n, C) as above, 
consider the following conditions: (Ml) Xk &U, k^l,. . . ,N; (M2) if rk> I, then h{z) is 
analytic in a neighborhood Uk 9 Afc, where Uk ^ U. Let Mh = {A G A/(n,C); conditions 
(Ml) and (M2) hold at the characteristic roots Ai,...,AAr of A}. We now define the 
primary matrix function h{A) associated with the scalar- valued stem function h{z). 



Definition 3.1. The primary matrix function h : A4h — > M{n, C) is defined as 



whe 



h{K) = Ph{J)p-^ ^ P 



h{Jx 



( HXk) 

h'iXk) 
h^'-'-'\Xk) 





h{Xk) 







p- 



\ 





h'{Xk) h{Xk) 



The following technical result is proved in Appendix A. The functions (t — w)^, r(D-|- 
/), eT'"8n(2:)i7tD/2 appearing in the result below arc all primary matrix functions. 

The same interpretation is also adopted throughout the rest of the paper, for example, 
with functions sin(7TD/2), cos(7tD/2) appearing in Theorem 3.2. (It should also be noted, 
in particular, that the definition of the matrix exponential based on a series is equivalent 
to that based on primary matrix functions.) 
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Proposition 3.1. Under (3.1) and condition (3.18) in Theorem 3.2, 

e'"=^((< - u)^ - i-u)^) du = V r^rP + /)eT^'S'^(")'"^/2. (3.17) 



Next, we construct time domain representations for OFBMs, which provides the main 
result in this section. Further comments about the resuh can be found after the proof. 

Theorem 3.2. Let {BH{t)}t£m be an OFBM with o.s.s. exponent H having the spectral 
representation (3.2) with A~ Ai +iA2, where Ai,A2 G A/(n,M). 

(i) Suppose that H Cz AI{n,R) has characteristic roots satisfying (3.1) and 

Re(/ife)^^, fc = l,...,n. (3.18) 
There are then Af-|_,Af_ e il/(n,R) such that 
{BH{t)}tm 

(3.19) 

= [jm - «)? - {~u)^W+ + {{t u)^ (-u)^)M^)B{du) 

where {_B(u)}„gK 'is a vector-valued process consisting of independent Brownian 
motions and such that EB(du)B{du)* = du. Moreover, the matrices M_ can 
he taken as 

M± = y|(^sin(^^) \{D + I)-'A,±cos(^^^ 'viD + iy^A^y (3.20) 

(ii) Suppose that H = (1/2)/. There then exist M,N G il/(n,IR) such that 



{BH{t)]teVL={ I ( (sign(<-u)-sign(-M))Af + log( )iV)B(du) 

(3.21) 

where {B(u)}u£m is as in (3.19). Moreover, the matrices M, N can be taken as 

M^^Ai, N = -^A2. (3.22) 

Proof, (i) Denote the process on the right-hand side of (3.19) by Xh- By using the 
Jordan decomposition of D, it is easy to show that Xh is well defined. It suffices to show 
that there are M± such that the covariance structure of Xh matches that of the OFBM 
Bh given by its spectral representation (3.2) with A = Ai +iA2. By using the Plancherel 
identity, note first that 

EXH{s)XHity 



\t~u\ 
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1 f (c'"^-l)(c-'*^-l) 



2n 



Meanwhile, for Bh, we have 

EBH{s)BH{ty^ / ^ ^(x7^AA*a;7^ +.tI^AA*x:^ )dx. (3.23) 

Thus, by using the relation e'^ = cos(0) + isin(0), Q £ M{n), it is sufhcient to find 
M± e M(n,R) such that 

AA* = ^T{D + /)(c"'"^/2a/+ + e'"^/2^/_) 

ZTT 



X (M;e'"^'/2 + Ar e-'"^'/2)r(L' + /)* 

^r(i:) + /) ^sin (^^^y) " ~ M* ) sin 



ni 

V 2 ; 



cosf ^) + + Ml) COS f'^'^* 



(3.24) 



i I cos I ^) + Af_)(M; - ir ) sin f 



- sin (^^) (A/+ - A/_)(A/; + Ar ) cos (^^^ ) ) r(i? + /)* . 

On the other hand, 

AA* = (Aivlt + A2A;) + i{A2Al - AiA*2). (3.25) 

By comparing (3.25) and (3.24), a natural way to proceed is to consider Af+ and A/_ as 
solutions of the system 

A, = -^TiD + I) sU^) (A/+ - Af_), 

(3.26) 

A2 = -^r{D + I)cos(^]{M++M^). 



2n " ' V 2 

By assumption (3.18), sin(^), cos(^) and r(_D + /) are invertible, and we obtain the 
solution given by (3.20). 

(ii) In this case, one can readily compute the inverse Fourier transform of the integrand 
in (3.2), that is (up to (271)-^), 



e'f^ - 1 



^ (1{2;>0}^ + l{a:<o}^) dx 
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f / cos((t — u)x) ~ cos((t ~ u)x) + i(sm((t — u)x) + sm(ux))\ , , -r, , 

= J{ — — ^^j(l{,>o}A+l{,<o}^)dx. 

As shown in Appendix B, this becomes 

-21og( )A2 + (sign(t - u) - sign(-u))7t^i. 



Then, by considering second moments and using Plancherel's identity, representation 
(3.21) holds with M = (27t)'1/27tAi and N ^ {2n)~^/^{-2)A2. It is weh defined because 
the integrand comes from the inverse Fourier transform of a square-integrable function 
and hence is also square-integrable. □ 

Remark 3. 3. Note that the invertibility of M or N in (3.19) is not a requirement for the 
process to be proper. A simple example would be that of a bivariate OFBM Bh whose 
time domain representation (3.19) has matrix parameters H = hi , h E (0, l)\{l/2}, 

1^0 0^' ^^0 1 

The two components of Bh are two independent (univariate) FBMs with exponent h. 
Thus, Bh is proper. 

Example 3.1. When (3.18) does not hold and H ^ (1/2)7, the general form of time 
domain representations can be quite intricate. For example, with 

D=i' 0\ / _^l/2 



1 oy V V 1 1/2 

the calculation of the inverse Fourier transform (up to {2n)^^) 



^{x^^A + xZ^A)dx (3.27) 

in Appendix B shows that Bh has the time domain representation 

{BH{t)}tm={ [ {fi{t,u)M + f2it,u)N)B{du)] , (3.28) 
UR J teR 



where M ^i, N 



sign(t — 7i) — sign(— u) 
(C + log \t — u|) sign(f — u) — (C + log |u|) sign(— u) sign(t — u) — sign(— u) 

/ lo6(ti) \ 



\og{^^] C+Uog{\t-u\\u\)) lo/^' 



2 / \ m J / 
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where C is Euler's constant. Note that, without taking the Fourier transform of (3.28), 
it is by no means obvious why its right-hand side has stationary increments and is o.s.s. 



4. Conditions for properness 

We now provide sufficient conditions for a process with spectral and time domain repre- 
sentations (3.2) and (3.19), respectively, to be proper and, thus, to be an OFBM. 



Proposition 4.1. Let {X(t)}t^R be a process with spectral domain representation (3.2), 
where the characteristic roots of H satisfy (3.1). If YLe{AA*) is a full rank matrix, then 
{X(t)}tgR is proper (i.e., it is an OFBM). 



Proof. We must show that 

tT 1 2 



Ex{t)x{ty 



IX 



{x^'^AA*x+^' +xZ'^AA*xZ'^'')dx, t^O, 

is a full rank matrix. For simplicity, let d/i(.T) — | ° .^7^ P d.T. Then 

EXit)Xit)* = [ x^^AA*x+"' d^iix) + [ x+^AA^x^^' d/i(a;) 

= 2 [ x-^Rc{AA*)x-^' dfi{x). 

The matrix J^x^^ Iic(AA*)x^^ d^{x) is Hermitian positive semidefinite. Moreover, for 
any vGC"\{0}, 

where the strict inequality follows from the fact that {v*x^^)Kc{AA*){x^^ v) > for 
all X > 0, the latter being a consequence of the invertibility of x^^ and the assumption 
that Re(AA*) has full rank. □ 



Based on Proposition 4.1, we can easily obtain conditions for properness based on 
time domain parameters. Consider a process {X{t)}ti:m with time domain representation 
(3.19), where the characteristic roots of H satisfy (3.1) and (3.18). If 

are full rank matrices, then {X{t)}tes. is proper (i.e., it is an OFBM). 



Remark 4-1 ■ Re(AA*) having full rank does not imply that A A* has full rank since 
1(^2^^ — A1A2) may have negative eigenvalues. Also, note that Ke{AA*) being a full 
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rank matrix is not a necessary condition for properness. For example, consider the process 
{X{t)}t£R with spectral representation (3.2), where 



AA* 

Then 
EX{t)X{t)* = 



1 2 

2 4 



H - 



hi 

/l2 



/ll,/l2€ (0,1). 



2 / |^|_2(hi-l/2) 2|2;|-((^i-l/2) + (?i2-l/2))- 

2|a;|-(('u-l/2) + (/i2-l/2)) 4|^,|-2(h2-l/2) 



<P''^C2(/ll)2 



where 



hi + /i2 



C2(/l)' 



4|t|2'^^C2(/l2)' 



7t 



/ir(2/i)sin(ft.7T) 

(see, e.g., [31], page 328). Therefore, det{EX{t)X{ty) ^0 if and only if 

^1 + ^2^ ^ 



C2{hirc2{h2r - a 



dx 



(4.1) 



which generally does not hold. 



5. Time-reversibility of OFBMs 

We shall provide here conditions for an OFBM to be time-reversible. Recall that a process 
X is said to be time-reversible if 

{X{t)}teM^{X{-t)}teR. (5.1) 

When X is a zero-mean multivariate Gaussian stationary process, (5.1) is equivalent to 

EX{s)X{ty ^EX{-s)X{~ty*, s,tGR, 

which, in turn, is equivalent to 

EX{s)X{t)* =EX{t)X{s)*, s,teR. 

The next proposition provides necessary and sufficient conditions for time-reversibility 
in the case of Gaussian processes with stationary increments. It is stated without proof 
since the latter is elementary. 
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Proposition 5.1. Let X be a Gaussian process with stationary increments and spectral 
representation 

c [ f e"^-l- 



{X(t)}teK- / —. Y{dx 



where Y{dx) is an orthogonal-increment random measure in C". The following state- 
ments are equivalent: 

(i) X is time-reversible; 

(ii) EY{dx)Y{dx)* = EY{-dx)Y{-dx)* ; 

(iii) EX{s)X{t)* = EX{t)X{s)* , s,t€R. 

The following result on time-reversibility of OFBMs is a direct consequence of Propo- 
sition 5.1. 

Theorem 5.1. Let {BH{t)}t£Ts. be an OFBM with exponent H and spectral representa- 
tion (3.2). LetA = Ai+iA2, w/iere Ai, A2 € M(n,R). Then Bh is time-reversible if and 
only if 



AA*^AA* or A2A\^ A^A*^. (5.2) 

Proof. From Proposition 5.1(ii), time- reversibility is equivalent to 

E{{x-^^A + x-_°^)B{dx)B{dx)*{A*xZj^' + 'AxZ^')) 

= E{{x-_°A + x-^°J.)B{-dx)B{-dx)* {A* xZ^' + A^a;;^* )) 



Since \x\ is invertible for a; > 0, this is equivalent to (5.2). □ 

Corollary 5.1. Let {-Bif(i)}teK be an OFBM with time domain representation given by 
(3.19) and exponent H satisfying (3.1) and (3.18). Then Bh is time-reversible if and 
only if 

cos ( ^ ) ( A/+ + ) (M; - M*_ ) sin 



2 

(5.3) 

sin [ — ) (A/+ - M_ ) (MT + A/* ) cos ( 



Proof. As in the proof of Theorem 3.2, under (3.18), the matrices sin(7tL)/2), cos(7tZ?/2) 
and T{D -\- 1) are invertible and, thus, by using (3.26), one can equivalently re-express 
condition (5.2) as (5.3). □ 



18 



G. Didier and V. Pipiras 



A consequence of Theorem 5.1 is that time-irreversible OFBMs can only emerge in the 
multivariate context since condition (5.2) is always satisfied in the univariate context. 
Another elementary consequence of Proposition 5.1 is the following result, which partially 
justifies the interest in time-reversibility in the case of OFBMs. 

Proposition 5.2. Let {BHit)}tm be an OFBM with H satisfying (3.1). If {BH{t)}tm 
is time-reversible, then its covariance structure is given by the function 



EBH{s)BH{tr = -(|t|^r(l, l)\tr + \sfT{\, l)\sf' 

-\t-sfT{i,i)\t^sr), 



(5.4) 



where r(l, 1) = EBh{1)Bh{1)* . Conversely, an OFBM with covariance function ( 5.4 ) is 
time-reversible. 



Proof. This follows from Proposition 5.1(iii) 



□ 



Remark 5.1. One should note that for a fixed exponent H, not every positive definite 
matrix r(l, 1) leads to a valid covariance function (5.4) for time-reversible OFBMs. 
In fact, fixr(l,l) = /, n = 2. We will show that, for an exponent of the form 



H = 



h 
1 h 



he [0,1), 



there does not exist a timc-revcrsiblc OFBM Bh such that EBh{1)Bh{1) 
From Theorem 5.1, 



EBHil)BHiiy 



where (s,j)i.j=i.2 := AA* e M{n,\ 



1 



\x\-''AA*\x\ 



-D* 



-2d 



IX 

Wc have 
1 



log(x) I) \Si2 S22 



Sll Si2 



-D* 



dx. 



r(i,i). 

(5.5) 



1 log(a:) 
1 



Sll Siilog(x) -I-S12 

Sll log(a;) 4- si2 sii(log(x))2 H- 2si2 log(x) + S22 



For notational simplicity, let 

We obtain 

EBH{l)BH{ir 



IX 



{\og{x)f\x\-'^'^ dx 



k = Q,l,2. 



siiro(rf) siiri(d) -|-si2ro(d) 

siiri(d) -f si2ra{d) siir2{d) + 2si2ri{d) 4- S22ro(d) 
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On the other hand, for any real symmetric matrix, the condition for it to have equal 
eigenvalues is that the discriminant of the characteristic polynomial is zero. In terms of 
£'i?//(l)i?fl'(l)*, this means that 

sii'^o('^) = siir2{d) + 2si2ri{d) + S22roid), siiri(d) + si2ro(d) = 0. 

Therefore, sn = si2 = S22 = 0, which contradicts the assumption that r(l, 1) = /. 

This issue is a problem, for instance, in the context of simulation methods that require 
knowledge of the covariancc function. For time- reversible OFBMs with diagonalizable H, 
one natural way to parameterize r(l,l) is by means of the formula (5.5) since, in this 
case, the former can be explicitly computed (see [16]). 

Finally, we provide a result (Proposition 5.3) characterizing time-reversibility of some 
OFBMs in terms of their symmetry group Gi (see Section 2.2). This result will be used 
several times in the next section. 

Proposition 5.3. Let {-B_ff(t)}tGR be an OFBM such that hi G £{Bh) for some h G 
(0,1). Then {-Bif(i)}teK is time-reversible if and only if Gi{Bh) is conjugate to 0{n). 

Proof. Regarding necessity, note that if such Bh is time-reversible, then, by Theorem 
5.1, its covariance function can be written as 

for some positive definite S € M{n) (note that if S is only positive semidefinite, then the 
process is not proper). 

For sufficiency, consider the covariancc function of the OFBM with exponent H ~ hi, 
h&{0,l), 

Define Bh = W-^Bh, where WO{n)W-^ = Gi{Bh) for a positive definite W. Then, for 
any O e 0{n), 

{OBH{t)}tm^{BH{t)}tm- 

By the uniqueness of the spectral distribution function, this implies that 0{W'^^ AA*W^^)0* = 
W-^AA*W-^, that is, 0{W-^AA*W-^) = {W-'^AA*W-^)0. Since O is any matrix in 
0(n), it follows that W~^AA*W~^ =cl,c€ C\{0} (for a proof of this technical result, 
see [8]). Thus, AA* = cW^ and c> 0. Hence, AA* = AA*. □ 
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6. The dichotomy principle 

We now take a closer look at the increments of an OFBM, which form a stationary 
process. 

Definition 6.1. Let {BH{t)}t&L be an OFBM. The increment process 

{YH{t))teT^{BH{t + l)-BH{t))teT, where T = Z or M, 
is called operator fractional Gaussian noise (OFGN). 

From Theorem 3.1, the spectral representation of OFGN in continuous time is 



{YH{t)}te^^\ I e"-^^(a;;^A + xI^A)i3(dx-) 



The spectral density of {YH{t)}t£VL is then 
|ei^_ 112 



/y^(x)= ' , ' (x-7"AA*a;7" + xZ"" AA* xZ"" ) , x€ 

since the cross terms are zero. 

In discrete time, analogously to the univariate expression, 



EYH{0)YH{n)* 

/n oo 
^inx fY^{x + 2nk)dx, neZ. 

-7T J. 



(6.3) 



fc=-c 

The spectral density of {YH{n)}n£i, is then 
9Yh{x) = 2(1 -cos(a;)) 

X E b + 27rfcp (^^ + ^^^)+''-^-^*(" + ^^^^+''' (^-4) 

k— — oo 

+ {x + 2nk)Z^AA*{x + 2nk)z'^'), xe [-7i,n]. 
The form (6.4) of the spectral density leads to the following result. 

Theorem 6.1. Let H be an exponent with (possibly repeated) characteristic roots hi, 
/ = 1, . . . , n, such that 

l/2<Re(/iO < 1, Z = l,...,n. (6.5) 

Let gvHi^) — {.9iff (2;)ij} be the spectral density (6.4) of OFGN in discrete time. Then, 
for fixed i,j , either: 
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(i) \gYH{x)ij \ ^ oo, asx^O; or 

(ii) gvaix)!] =0 on [-7T,7t]. 

Proof. Let di and N be as in (2.2) and (2.3), and take x > 0. By assumption (6.5), 
< Re((i;) < 1/2. For a given z > 0, if we take —D in Jordan canonical form PJP^^, we 
obtain that 

z"^ = Pdiag(z'^-''i , . . . , z'^-''«)p-\ 

where J_d, is a Jordan block in J, I = 1, . . . , N < n. Without loss of generality, for A: > 0, 
each term of the summation (6.4) involves the matrix expression 

Pdiag((a; + 2nk)i-'" , . . . , (x + 2nk)i-'''^)p-^A 

(6.6) 

X A*{P*)-^ diag((a; + 27Tfc)+*''i , . . . , (x + 27tfc)+"''" )P*. 

Denote the entries of the matrix-valued function (6.6) by h{x + 2nk)ij, i,j = l,...,n. As 
shown in Appendix D, h(x + 2nk)ij is a linear combination (with complex coefficients) 
of terms of the form 

p{x + 2nk){x + 2'nky^'q{x + 2nk){x + 27Tfc) /, m = 1, . . . , A^, 

where p{x),q{x) are polynomials (with complex coefficients) in log(x). Thus, 



sup 



,X + 27tfcp 

fe=l ' ' 



<oo, i,j = l,...,n, (6.7) 




and, therefore. 



hm 2(1 - cos(x)) yj^h{x),j + }^ J—^r^^H^ + '^^^h 

= lim 2(1 — cos(x)) I -r-T^h{x)ij 
x^o \ \xY 

On the other hand, since Re(d;) > for / = 1, . . . , n, h{x)ij diverges as the power function 
(times some p{x)q{x)) as x — > unless it is identically zero for all x (in particular, for 
X + 2nk). Thus, by (6.7) and the fact that ^^^"^"'^'■^^^ — 1 as x — 0, the claim follows. □ 

In the univariate context, the range (1/2, 1) for H is commonly known as that of long- 
range dependence (LRD). In the multivariate context, characteristic roots of H with real 
parts between 1/2 and 1 have the potential to generate divergence of the spectrum at 
zero. Theorem 6.1 thus states that if OFGN is long-range dependent in the sense of (6.5), 
then the cross correlation between any two components is characterized by the following 
dichotomy: 
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• it cither has a divergent spectrum at zero, a characteristic usuahy associated with 
LRD; or 

• it is identically equal to zero. 

Obtaining a similar dichotomy principle for a larger range of characteristic roots than 
that in (6.5) is much more delicate. The following two examples illustrate some of 
the potential difficulties. Example 6.1 shows that if one of the characteristic roots of 
D = H — (1/2)/ is 0, then the dichotomy may not hold. Example 6.2 shows that cer- 
tain cancellations may occur in the cross spectrum if the characteristic roots of D have 
opposite signs. 



Example 6.1. If, for instance, D = P dia.g{d,0)P ^, where < c? < 1/2, 

P 

and A := P, then, as a; — > 0, 



1 V2/2\ 
v^/2 J 



, . /a;-2'i + i/2 1/2 

where '-^ indicates entrywise asymptotic equivalence. As a consequence, if one of the 
components of OFBM behaves like Brownian motion, then this may create cross short- 
range dependence among the components. This example is a direct consequence of a 
more general operator parameter D whose eigenspaces are not the canonical axes. If we 
take, instead, D = diag((i, 0), whose eigenspaces are the canonical axes, each term of the 
summation (6.4) has the form 

diag((:r + 27rfc)±^ 0) AA* diag((x + 27tfc)±^ 0) = + 27tfc)±2'i o 

where 

(s«j)«j=i,2:=^^* (6.8) 

and thus the dichotomy holds. 

Example 6.2. Consider Ae GL{n,R) and D = diag(d, -d), where de (0,1/2). Using 
the notation (6.8), we have 



gvH (x) - Pdia.g{x-'^,x'^)P*AA*Pdiag{x-'^, x'^)P* 



SllX ^'^ Si2 
Sl2 5222:^'' 



as a; —7- 0. Here, the multivariate differencing effects of the operator D cancel out in the 
cross-entries. 



From a practical perspective. Theorem 6.1 raises the question of whether the class 
of OFGNs is flexible enough to capture multivariate LRD structures. This, and related 
issues regarding multivariate discrete time series, will be explored in future work. 
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7. Operator Brownian motions 

In this section, we shall adopt the following definition of multivariate Brownian motion 
and establish some of its properties. 

Definition 7.1. The proper process {BH{t)}te9. '"^ operator Brownian motion (OBM) 
if it is a Gaussian o.s.s. process which has stationary and independent increments and 
satisfies Bh{0) ~ a.s. 

In place of the condition Bh{Q) = a.s., we can assume that the characteristic roots 
of the o.s.s. exponent H have positive real parts, which implies the former condition. 
Another important way to motivate Definition 7.1 is as follows: since {i?jy(i)}tgR is L^- 
continuous (see the beginning of the proof of Theorem 3.1), by Hudson and Mason [18], 
Theorem 4, and Hudson and Mason [19], Theorem 7, our Definition 7.1 implies that 
(1/2)/ can always be taken as an exponent of OBM. 

The next proposition and example show that an OFBM Bh with (1/2)/ e £{Bh) is 
not necessarily an OBM. This stands in contrast with the univariate case. 

Proposition 7.1. Let {BH{t)}t^R be an OFBM with exponent // = (1/2)/. Consider 
its time domain representation (3.21) with parameters M and N , or its spectral domain 
representation (3.2) with A = Ai+ 1A2. Then {/?ir(t)}teK is an OBM if and only if the 
following two equivalent conditions hold: 

(i) MN* = NM* ; 

(ii) A2Al=AiAl. 

Proof. Since (1/2)/ G £{Bh), it follows that Bh{0) = a.s. Therefore, we only have to 
establish that the increments are independent if and only if (i) holds. Demonstrating the 
equivalence between (i) and (ii) is straightforward by using the relation (3.22). 
Write the time domain representation (3.21) of Bh as 

/ (2(l{t_„>o} - l{-„>o})M + (log \t^u\~ log I - u\)N)B{du). (7.1) 
Js. 

Take si < < S2 < ^2- For the increments of the process Bh, we have 

E{BH{ti) - BH{si)){BH{t2) - Bh{s2)T 

= / (AA{,,<^^t,^l{s,<u<t,}MM*+\ogf^^\ogf^^NN* (7.2) 
+ 2 • l{s,<u<t,} log i^i^A/iV* + 2 • l{s,<u<t,] log i^l^iVA/*) An. 

\S2-U\ \S1-U\ J 



From the univariate time domain representation of Brownian motion, the first two of the 
four terms in (7.2) have zero integral. Define 'p{u) = u(log(M) — 1). The right-hand side 
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of the expression (7.2) then becomes 



{ip{t2 - Si) - ip{t2 - h) - ip{s2 - Si) + Lp{s2 - ti))MN* 

+ {ip{t2 - h) - ip{s2 - h) - ip{t2 - Si) + (p(s2 - Si))NM*, 



which is identically zero if and only if M N* = NM* . 



□ 



Example 7. 1 . Consider a bivariate process X defined by the expression (7.1). Set M = / 

and let iV = L e so(2)\{0}. Then MN* ~NM* ^ (from which the cross terms in 
expression (7.2) cancel out when si = S2 = and ti=t2~ t) and 



which is a full rank matrix for t ^ Q (to obtain the constant tt^, one can use, e.g., 
Proposition 9.2 in [33] and Proposition 5.1 in [32]). Hence, X is proper. This gives an 
example of an OFBM for which (1/2)/ € £{X) but which is not an OEM. Moreover, it 
is an example of an OFBM with an exponent of the form /i/, h € (0, 1), but which is not 
time-reversible and for which Gi = 0(2) docs not hold by Proposition 5.3 (in contrast, 
by Hudson and Mason [19], Theorem 6, Gi{X) = 0{n) implies that hi g S{X) for some 



The following result is a direct consequence of Theorem 5.1 and Proposition 7.1. It 
shows that in the class of OFBMs with exponent H = (1/2)/, time-reversibility is equiv- 
alent to independence of increments. 

Corollary 7.1. Let {/3//(t)}tgM be a time-reversible OFBM with exponent H = (1/2)/. 
Then {BH{t)}teR *s an OEM. Gonversely, let {S_f/(i)}tgR be an OEM. It is then time- 
reversible (and has exponent H = [1/2)1). 

Remark 7.1. Note that, as a consequence of Proposition 5.3, time-reversibility may be 
replaced in Corollary 7.1 by the condition that Gi is conjugate to 0{n). In other words, 
an OBM is elliptically symmetric. 

We conclude by providing a spectral representation for OBM. 

Proposition 7.2. Let {BH{t)}t(iR be an OEM. Then 



for some positive definite operator W , where {i?(i)}tgR is a vector of independent stan- 
dard EMs. 

Proof. Consider the spectral domain representation of Bh with parameter A = Ai+ \A2 . 
Set W := {AiAl -1-^2^2)^/^. The result follows from Proposition 7.1(ii), and relations 



EX{t)X{t)* = A\t\I + TT?\t\L{-L) = 1*1(4/ - n^L^) 



h). 




(7.3) 



(3.23) and (3.25). 



□ 
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Appendix A: Fourier transforms of OFBM kernels 



In this appendix, the goal is to prove Proposition 3.1. First, we state a condensed version 
of Horn and Johnson [17], Theorems 6.2.9 and 6.2.10, pages 412-416, which will be useful 
in the subsequent derivations. We shall use the notation introduced before Definition 3.1. 

Theorem A. 1. Letf,g:U^C be two stem functions and let A4 fg ~ A4 fOMg. Then: 

(i) the primary matrix function f :Aif ^ M{n,C) is well defined in the sense that 
the value of /(A), A G A^/, is independent of the particular Jordan canonical 
form (i.e., block permutation) used to represent it; 

(ii) /(A) = g(A) zf and only if /(^^ (Afe) = g^^' (A,.) /or j = 0, 1, . . . , r,. - 1, A; = 1, . . . , TV 
and K & M fg; 

(iii) for q{z) := f{z)g{z), we have q{A) = /(A)5(A) = 5(A)/(A) for A e Mfg; 

(iv) for s{z) := f{z) + g{z), we have s(A) = /(A) + g{h) for Ae Mfg. 

Throughout this section, we assume (3.1) and (3.18). Denote by F the Fourier trans- 
form operator. For c? G C such that 

Re(rf)e(-l/2,l/2)\{0}, (A.l) 

define 

f±{t,u,d)^{t~u)i^{-u)i 

and 

pita; _ 1 

/i±(i,x,d) = |a;|-'^r(d+l)eT'^'sn(.)i7td/2^ 

It is well known that 

T{f±{t,-,d)){x)^h±{t,x,d) (A.2) 

when de (-1/2, l/2)\{0} (see, e.g., [29], page 175). One can show that (A.2) also holds 
under (A.l) (see Remark A.l). 

For the purpose of calculating Fourier transforms of primary matrix functions associ- 
ated with the stem functions f± and h± , we will need to consider derivatives of the latter 
with respect to d. Note that, for fixed x, the functions T{d-{- 1), cTsiga(x)ind/2 ^^^^ \x\~'^ 
are holomorphie on the domain —\< Re(d) < ^. Thus, so are the functions h±{t,x,d). 
Note that, for fixed t and u, since (t — u)^ and (— u)^ are holomorphie on the domain 
-1/2 < Rc{d) < 1/2, then so are f±{t,u,d). 

As a consequence, by Theorem A.l(i)-(iv), 

/l±(t,X,i:>) = i|x|--°r(D + /)eT'''sn(^)mD/2 

ix 

and 



/±(t,u,I?) = (t-u)±-(-")±. 
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We now need to show that 



■^(/± (t, u, D)){x) = h± {t, X, D), (A.3) 
where T is the entrywise Fourier transform operator. 



Proof of Proposition 3.1. We wih break up the proof into three cases: 

Case 1: —1/2 < Kc{dk) < 0, k = 1, . . . , N . We wiU develop the calculations for h^^, 
which can be easily adapted to h^. By Theorem A.l(ii), in the case of h^^, (A.3) is 
equivalent to 

— h+{t,x,d) = -g^ jj^-f+{t,u,d)du = jj^-^f+{t,u,d)du (A.4) 

at d = dfc, for fc = 1, . . . , A^, j = 0, 1, . . . , Tfc — 1. Consider the domain A{d,d) {dGC:d< 
Re(d) < d)}, where —1/2 < d < d < 0, which is open and convex. Consider j = 1, that 
is, the first derivative. Fix d* G A{d,d) and take a sequence {dmjmeN C A{d,d) such 
that dm — > d* . For each m, by the mean value theorem for holomorphic functions ([13], 
Theorem 2.2), there exist constants 5i{m). where 

5i{m) = am.idm + (1 - am,i)d* , arn,i S (0,1), i = 1,2, 



such that 

f+{t,u,djn) - f+{t,u,d*) 
d„, — d* 



Rc 



f+(t,U,djn) - f+{t,u,d*) 



i Im 



d,n - d* 
f+{t,u,dm) - f+{t,u,d*) 



d* 



d 



dd' 



Rc — /+(t,u,5i(m)) +ilm — /+(t,u,52(m)) 



d 



dd' 



We will now obtain an integrable function that majorizes \-^fj^{t,-,d)\ for all c? € 
A(d, d). Assume, without loss of generality, that i > 0, and take d € A(d, d) and 5 > 
such that —1/2 < Re(d) — 5 and Re(d) + <S < 0. From the continuity of -ggf+{t,u^d) for 
< u < t, there exist constants Ki and rji such that 



<\\og{t-u)+\\{t-u)'i\ 
<Kr^,,t-rn]{n) + \{t-u)l'^^^''\ 



(t-nut) 



(A.5) 



u>0. 



Also, there exists a constant K2 such that 

I log(^ - u)+{t - uf^ ~ \og{-u)+{-uf^\ <K2 + |(-u)f -1 < < 0. (A.6) 
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One can show that there exist constants A'3 and 7^2 < ~1 such that 

\\og{t-u) + {t-u)i^\og{^u)+{~u)i\<K3{-u)f~'^+'-\ u<m- (A.7) 

From (A. 5), (A. 6) and (A.7), and from the fact that u, d) is bounded on 772 < u < 

— 1 uniformly in d € A(d, d), we conclude that the ratio ■^+(*''''^^)~-^+'^*^'^'^ ) bounded 

by a function in L^(R). Thus, by the dominated convergence theorem (for C- valued 
functions), we have 



. f+{t,u,dra) - f+(t,u,d*) 
dm - d* 



du - 



■ d 

e""— /+(t,u,(i*)du, m^oo. 
od 



We can always assume that t 7^ and the case of t < can be dealt with in a similar 
fashion. The extension of the above argument for derivatives of higher order j poses no 
additional technical difficulties. This establishes (A. 4). 

Case 2: < Yic{dk) < 1/2, k = 1, . . . , N . In this range, the upper bound in (A.7) is not 
in L^(K) so we need a slightly different procedure. Since h±{t, •, d) € L^(R). we can apply 
on both sides of equation (A. 2) and obtain 

f±{t,U,d)^T-\h±it,;d)){u). 

Therefore, it suffices to show that 

f±it,u,D)=T-\h±{t,x,D)), 



(A., 



where ^ is the entrywise inverse Fourier transform. 
Note that expression (A. 8) is equivalent to 



gJ-f±{t,X,d)^ — 



'h±(t^ X, d) dx -■ 



-^h±{t,x,d) dx 



at d = dk, for fc = 1, . . . , A^, j = 0, 1, . . . , — 1. To show this, one may proceed as in the 
case of — 1/2 < Re((i) < 0. The existence of an upper bound in L^(R) is ensured by the 
fact that 



ddJ 



1 



< 



< 







\ ix 











Atx 



log|x|P>|-^'=('') 

iog|xii-'>ri^^®i{o<i.i<i} 

l0g|x|P>|-«'°(^^)l{l<|,|<^} 



(A.9) 



which is integrable for all d € A(d, d), < d < d < 1/2. 
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General case: As a consequence of (A. 3), 

T{f±it,;J)){x) = h±{t,x,J) (A.IO) 

holds, where J is a matrix in Jordan canonical form with characteristic roots satisfying 
(3.1) and (3.18). Now pre- and post-multiply equation (A.IO) by P and P^^, respectively. 
Since 

Pr(- J)P-1 = r{-D), Pei^^/2){J+I)p-l ^ ^(in/2)(D+I) ^ 

it follows that 

Pf±{t,u,J)p-^ = f±{t,u,D), Ph±{t,u,J)p-^ = h±{t,u,D), 
from which we obtain equation (A. 3). □ 

Remark A.l. A common way to prove that (A. 2) also holds for d£C satisfying (A.l) 
is by analytic continuation. In particular, this requires the ability to differentiate under 
the integral sign in the Fourier transform. The latter could be achieved by following the 
argument in the proof of Proposition 3.1. 



Appendix B: Some useful integrals 

In this appendix, we calculate the inverse Fourier transforms used in the proof of Theorem 
3.2(ii) and Example 3.1. We shall use several formulas from [14]: 

f ^ sm(ax) , f ^ sin(ax) , n . 

/ 1{:e<o} Ax^ / l{^>o} da; = -sign(a) (page 423), (B.l) 

cos(aa;) — cos(6a;) , , 161 , , ,^ 
1{.>0} ^da; = log^ (page 447) (B.2) 

( , , f cos(aa;) — cos(6a;) , , l&l \ ,^ „n 

therefore, / l^^^^^}^ ^dx^-log^ , (B.3) 

' dx 7T 

log(2;)sin(aa;)— = (C + log(a)), a>0 (page 594) (B.4) 

X 2 



(B.5) 



Ax f dx 
therefore, with a e M, / log(a;) sin(aa;) — = / log(a;_) sin(aa;) — 

^ J — CO 

= -1^{C + \og\a\)s\gn{a) 



I log(a::)(cos(aa::) — cos(6a;))— = logf y I f C + i log(a6) I , a,6>0 
io X 2 ) 



(page 594), 
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where C is Euler's constant 



^therefore, J log(a::_)(cos(aa;) — cos(6a;))— = — log [ -j^ ] [ C + log(|a5|) 



(B.7) 

Using these formulas, wc obtain that, for a; > 0, 



gita: _ ^ 
^-iux — , l{^>o}dx 



— ^cos((t — u)x) — cos(ua;) + i(sin((t — u)x) + sin(wx))^ l{a;>o} dx 



1 , / |u| \ 7T . 

= T log ( i^Z"^ ) + 2 (^^S"^*-* - u) - sign(-M)). 



Similarly, for x < 0, 



e — 1 



— ^cos((t — u)x) — cos{ux) + i(sin((t — u)x) + sin(ita;))^ l{a;<o} dx 



1 , / lul \ 7T . 
-T log[ ^ I + ^(sign(t - u) - sign(-u)). 



Therefore, 



/ e-™"^ (l{,>o}^ + l{,<o}^)dx 

= (sign(i -u)- sign(~u))i Re(A) + log [t^A ^ 

which is the formula used in the proof of Theorem 3.2(ii). 

We now turn to the calculations of the inverse Fourier transform (3.27) in Example 
3.1. Note that 

We only need to calculate the inverse Fourier transform of the log term on the lower 
off-diagonal. 

For X > 0, using the formulas above, 

r Qit^ _ \ 

/ 0-'"=^^ log(x+)l|^>o}dx 
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1 



(cos((i — u)x) — cos{ux) + i(sin((t — u)x) + sin(ua;))) log(x+)l{2.>o} dx 



- ^((C + log \t - u\)sign{t -u)-{C + log |u|) sign(-M)). 



Similarly, for x <0 



f e 5-log(a;_)l{^<o}dx 



— (cos((i — u)x) — cos(ua;) + i(sin((t — u)x) + sin(ua;))) log(x_)l{a;<o} da; 



IX 

= -\^og(^-^){c+\^ogi\t-u\\u\) 

- ^((C + log \t - u\)sign{t -u)-{C + log \u\) sign(-u)). 

Therefore, for a G C, 

, e'*^ - 1 



IX 



(- log(a-+)l{3.>o}a - log(a:_ )l{^<o}a) dx 



= ((C + log|t - u\)sign{t -u)~iC + log|u|) sign(-w))i Re(a) (B.9) 

By combining (B.9) and (B.8), one obtains the time domain kernels on the right-hand 
side of (3.28). 

Appendix C: Nonexistence of OFBM for certain 
exponents 

Proposition C.l is mentioned in Remark 3.2. 

Proposition C.l. There does not exist an OFBM with exponent 



Proof. Assume that such an OFBM exists. For notational simplicity, denote the process 
by X, and its entrywise processes by Xi and X2. We will show that X is not a proper 
process. 
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Note that, for c > 0, from the matrix expression for and o.s.s 

c 



Xiict) 
X^ict) 



ten 



cXiit) 
clog(c)Xi(t)+cX2(t), 
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(C.l) 



In particular, this imphes that Xi is FBM with Hurst exponent 1. Thus, Xi{t) = tZ a.s., 
where Z is a Gaussian random variable (e.g., [33]). By plugging this into (C.l), we obtain 



ctZ 

c\og{c)tZ + cX2{t) 



(C.2) 



In particular, by taking c = t and t = 1, (C.2) implies that 

EX2{t)Z = E{t\og{t)Z + tX2{l))Z ^tlog{t)EZ^ + tEX2{l)Z. 

Thus, 

E{X2{t + h)- X2{h)){X^{l + h)- Xi{h)) 
= EX2{t + h)Z - EX2{h)Z 

^{{t + h) log(t + h) - h\og{h))EZ^ + tEX2{l)Z. 



(C.3) 



By stationarity of the increments, the expression (C.3) does not depend on /i, which is 
possible only when 

EZ'^ = 0. 

As a consequence, {X2{ct)}tem. = {cX2(t)}ieE- Thus, 



Xi{t) 





tY 



a.s. 



where y is a Gaussian random variable. In particular, X is not proper. 



□ 



Appendix D: The exponential of a matrix in 
Jordan canonical form 

Initially, let J\ £ M{n,C) be a Jordan block of size n\, whose expression is 



Jx 



/A 
1 A 
1 A 








(D.l) 



\0 ... 1 A/ 
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We have 



.Jx - 



(log z)z-^ 
2! 

(logz)"^-i 



V" 



1)1 





{log z)z'^ 

(logz)"^-2 

K-2)! ' 










(logz)z'^ z'^ 



(D.2) 



The expression for z'^, where J is, more generahy, a matrix in Jordan canonical form 
(i.e., whose diagonal is made up of Jordan blocks), follows immediately. In particular, 
the series-based notion of the matrix exponential is consistent with the primary matrix 
function-based notion of the matrix exponential. 
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